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Abstract 

We illustrate how one can easily derive kinetic energy operators for polyatomic molecules using 
polyspherical coordinates with very general choices for z-axis embeddings and angles used to specify 
relative orientations of internal vectors. Computer algebra is not required. 
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I. INTRODUCTION 


Mladenovic[l, 2] and Gatti et a/. [3-5] have presented derivations of kinetic energy op- 
erators for treating rotations and vibrations of polyatomic molecules using polyspherical 
coordinates. In these coordinates, the N nuclear position vectors are transformed to A — 1 
internal vectors and the nuclear center of mass vector, and the internal vectors are pa- 
rameterized by spherical polar coordinates: i = - 1- Then a body-fixed 

coordinate system is introduced, defined by aligning the body-fixed 2 axis along the fiist 
vector, and placing the second vector in the body-fixed xz plane. The remaining angles 
are either referenced with respect to the body-fixed z axis, or other vectors. [6] 

These coordinates are very useful, but are clearly not the most general choices. One choice 
of coordinates that is not treated is well known to be useful for triatomics. Sutcliffe and 
Tennyson [7] give the kinetic energy operator for triatomics using polyspherical coordinates 
in a very general form, and one finds that their expression having the body-fixed z axis lying 
between the two vectors to be very useful. For example, taking the body-fixed 2 axis to 
bisect the HOH angle in water gives quite good separation between rotation and vibration. 
One can see that this type of embedding might also be useful for systems of more than three 
atoms. Furthermore, for tetraatomics, it can be useful to have one of the internal angles 
giving the relative position of two light atoms, while another vector defines the body -fixed 
2 axis, e.g. having the HCH angle in formaldahyde a coordinate while the heavy CO group 
defines the body-fixed 2 axis. In addition, it may be desirable in this example to define 
the orientation of HCH as a unit vector relative to CO. While the first is possible using 
Mladenovic’s[6] local axis coordinates, the second is not. Ir. this work we will present a 
method for the derivation of kinetic energy operators for coordinates that encompass all of 
these choices. 

The method avoids complex intermediate triganometric expresions by factorizing out the 
dependence on Euler angles, which is an alternative to the scheme of Lukka,[8] and allows 
one to build up the kinetic energy operator for complex systems from fragements where the 
kinetic energy operator is already known. The treatment of non-orthogonal vectors is also 
straightforword as is the introduction of Eckart rotations. We only consider three, four, and 
five atom system in the present work, but the techniques described can be easily used to 
obtain the kinetic energy operator for more complex systems. 
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II. INTERNAL COORDINATES AND GENERAL PLAN 


Let X denote the 3 x N matrix of nuclear coordinates, and x the 3 x N matrix of internal 
vectors, with the last vector giving the position of the nuclear center of mass. Then these 
vectors are related by 

X = xM, (1) 


with M an N x N matrix of mass factors. We make two restrictions on M: first of all its 
inverse, denoted M, must exist, and secondly, there is no coupling of the other vectors to 
the nuclear center of mass vector. For a more in-depth discussion of these mass factois, see 
the work of Aquilanti and Cavalli,[9], Bramley and Carrington, [10] and Sehwenke.[ll] Then 
since the exact kinetic energy operator for the nuclei is 
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where m a is the mass of nuclei a , and i runs over x, y , and 2. we find that 
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where the reduced mass factors are given by 
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One can choose M so that — — is non-zero only when /? — /?, but that is not a requirement 

Pap* 

of our theory. So far we have done nothing new. 

We next re-write the kinetic energy operator in the symmetric form 
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where f means complex conjugate acts to the left. Since the center of mass decouples from 
the rest of the internal coordinates, in the sum we will take /? and /?' to run from 1 to A - 1. 
We then seek expressions of the form 
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with A a rotation matrix parameterized by the Euler angles a LB , ,8 LB , and 7 LS , which 
take the laboratory frame to the body- frame, g is a matrix depending only on vibrational 
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coordinates, and the O t are operators which include derivatives with respect to vibrational 
coordinates as well as total angular momentum components. Then, when we substitute this 
into Eq. 5, we obtain 

T = £d! GwOr, ( 7 ) 

iv 

where 

Gw = h 2 / 2 ^2 gjpi9ji3’i'/t*w- 
W 

Thus once we know g, we can very easily form G and hence the kinetic energv operator. 

To solve for ro-vibrational wave functions, we use the variational principle and analytic 
basis functions. Thus we can make explicit use of the symmetric form of the kinetic energy 
operator. For users of grid methods, it is required that both derivatives act to the right. 
However it is an easy excercise in differential calculus to convert the symmetric expressions 
to the more customary form. When doing this, one must bear in mind that in matrix 
elements of the kinetic energy operator of Eq. 7, none of the operators act on the Jacobian 
determinant. However, when one integrates by parts to obtain the expression with both 
derivatives acting to the right, the jacobian determinant is differentiated. 


III. ONE VECTOR 


Consider a single vector. Its spherical polar form in the space fixed axis system is 

/ R sin 9 cos 4> \ 
x = Rsin 9sm(f> 

\ R cos 6 } 

The associated angular mometum vector is given by 
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We will define a body-fixed z axis to be along f, thus to rotate vectors to this frame of 
reference, we will rotate by the rotation matrix A(<p, 9, 0) (see Edmonds[12] and Ref. 11 for 
rotation conventions), i.e. 
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and 
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L bf = A((f>,0,O )L = ih 

V o ; 

Now L bf is a rather particular angular mometum vector in that it does satisfy customary 
commutation relations, normal or anomolous, so one must exercise care in its use. One can 
easily show that 
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with gd) given in Table I, and O; made up of L b J , and L b J . This expression can be 
substituted into Eq. 5 to yield the kinetic energy operator: 
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Substituting Eq. 12 into this yields the well known kinetic energy for a diatomic written in 
symmetric form. 


IV. TWO VECTORS 

In this case Eq. 13 becomes 

( 15 ) 

9x v 7 b 

where O h comprises L b x [\ L b y { 1 , and O h comprises L b x { 2 , and L b y { 2 , where the 
superscripts on the angular momentum operators indicates they are expressed as in Eq. 12 
in the axis system particular to each vector. If our internal coordinates are orthogonal, 
i.e. 1/7 l Di3' = 0 when /? (3 1 , then we can substitute this into Eq. 5 and simply obtain 

the kinetic energy operator for two uncoupled diatomics. This is of course only of passing 
interest. To take advantage of the great simplification of total angular momentum coupling 
or to treat non-orthogonal coordinates, we need to remove the (3 dependence on the Euler 
angles in Eq. 15. When this is done, when we insert the expression for the derivative into 
Eq. 5, the rotation matrices vanish because of orthogonality. To do this we need to define a 
body-fixed coordinate system. This is done by a sequence of three rotations. We fiist rotate 
by A(g>i,0i,O), which takes the first vector to the body-fixed z axis and the second vector to 
the angles 9 2 and <p 2 - The second rotation is by A(0, 0, <p 2 ), which takes the second vector to 
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the x-z plane. The final rotation is A(O,a0 2 ,O), where a is a parameter running from zero 
to one. [7] When a = 0, the body-fixed 2 axis is along the first vector, while when a = 1, 
the body-fixed 2 axis is along the second vector, and when a = 1/2, the body-fixed 2 axis 
bisects the two vectors. 

We begin by considering only the first two rotations. We have 

A = £.4. lWl ,0,,O)£^>O,, (16) 

<*r x l j /, 

= E E °> &) E ( 17 ) 
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where 0 (l is the same as O tl except that the angular momentum operators have been rotated 
to the body frame of the triatomic. These will be denoted and Ly j . F 01 the second 
vector we have 
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To complete the transformation, we need to express the angular momentum operators 

in terms of -2- and the total angular momentum operators. We start by eliminating L^{ 
d&2 

and Lyf by introducing the total angular momentum vector J. We have to be careful here. 
In the space fixed frame, J = L\ + Z 2 . however we have L 2 in the frame BF obtained 
by rotating by A(<pi, 9\, ^ 2 )- If we rotate the space frame J by A(0i, 0i, <^ 2)5 obtain 
an angular momentum vector satisfying the anomolous commutation relations, thus it is 
customary to reverse the signs of the components, and we will do so here as well. Thus we 


have 
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Next we need Lff . Now L% 1 is given by Eq. 10 using <j> 2 and $ 2 as the independent variables, 
so we rotate L b 2 x bv A(0, 0, 0 2 ), which yields L^ 2 — f^cot an d ^y 2 = Finally 

we note that 

J, = ift-S-, (23) 
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and we use this to eliminate the derivative with respect to fa In Oi 2 , we need L i2 ~ ■ This is 
given by Eq. 12 using (f ) 2 and 0 2 as the independent variables. 

In summary" we now have 
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dxi 
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for fj =1 and 2, where is given by Table II with a — 0, and Oi comprizes ^ , a/? 2 ’ 
— J BF J BF . and J BF , where y = 0 2 is the angle between :he two vectors. We can then 
substitute this into Kcj. 8 ? and obtain the well known expression for the triutoinic kinetic 
energy. [13] 

For the most part, the introduction of the third rotation is straightforward. We insert a 
resolution of the identity so that Eq. 24 becomes 


Ajj(fa,0i,fa) A m j(0, o>0‘2 , 0)A mfc (0, ii02, 0) 9k0iOi 
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(25) 

(26) 

1 * ' 
with a LB , (3 LB , and q XB the new Euler angles relating the space fixed frame to the new 
body-fixed frame denoted BF. It is straightforward to multiply the matrices A and g (2) to 
get a new effective g^ 2 \ but more care needs to be taken to re..ate J in body frame BF to J 
in body frame BF. Formally we are changing variables from tq , 0\ , fa. and 0 2 to a LB , fl LB , 
7 iS , and y, where 

A(a LB , 0 LB , i lb ) = A (fa, On fa) A(0, a0 2 , 0), (27) 

and 

X = h. (28) 

Using the chain rule, we have 


d da LB d dp LB d , dp B d , d X LB d 


du> du da LB+ dio d/3 LB ^ du dj LB ^ dv d X LB ' 
oj = fa, $n fa, and d 2 . Now the only nonzero ^7- occurs for */ = 0 2: thus the equations for 
u 0 2 are just the equations describing a rotation by the constant a0 2 , thus 
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can be used to eliminate the BF components of the total angular momentum in terms of 
the BF components. This amounts to multiplying the final three columns of the result so 


far from the right by A7 (0. g$2j 0)- The final equation involving $2 is most easily analyzed 
considering the instantenous case where 4>i = = 4>2 = 0. Then we have 
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T- = -aiJ** + 
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dd 2 y dx 

i.e. the fifth column of the result so far needs to be aumented bv -ai times the third column 
Combining all of this together, we have 


d 

dxm 


'EM**™’! 31 ' 8 ’ 



(32) 


with g( 2 > given by Table II and O t comprizing ggp J** Jy T , and Jf*. Substituting 

this expression into Eq. 5 W'e obtain the expression of Sutcliffe and Tennyson]/] for theii 
z-axis embedding coordinates. 

It should be noted that g (2) in Table II is overly complex for triatomics. This is because 
we could have factored out the rotation matrix left multiplying g (2) , and this rotation matrix 
would collapse down to the unit matrix in Eq. 8. However this rotation matrix is required 
for the next sections, so we retain it here as well. 

To obtain the x-axis embedding of Sutcliffe and Tennyson, we simply insert a rotation 
that exchanges the x and 2 axes. This changes g (2) in only to a minor extent. We can also 
exchange the 2 and y axis, if desired. [14] 

Now consider the possibility that a depends on R\, R 2 i and x- The only changes are Eq. 


31 is replaced by 
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where the arrow means “replaced by”, and these changes modify g (2) in a straightforward 
manner. That is, it is not difficult to choose a so that the Eckart conditions[15] separating 
vibration and rotation are satisfied. The only complication is the kinetic energy operator 
matrix elements are no longer obtainable analytically, but it should be possible to obtain 
accurate results numerically with very little additional cost compared to calculations keep- 
ing a fixed. There is, of course, still the problem of singularities arising when the Eckart 
conditions lead to ill posed equations at linear configurations. 
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V. THREE VECTORS 


When we have three vectors, we use the results from the single vector case and the two 
vector case. We still use the notation for the total angular momentum in the final 
body-fixed frame, so we will change our notation for the two vector result. For the two 
vector sub-part, we will call the angular momentum j i 2 , the body-fixed frame defined for it 
bf 12, and the Euler angles og2, 0i2, and 7i2- Then w^e have 
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dx 


'i /3 
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for 0 = 1 and 2, with 0 (l2 comprizing ^ j b J 12 , J h y !u , and j b J n , and 
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with O h comprizing , and L b y { 3 . 

Next we have two choices. We can make the over-all Euler angles 0 3, 03 , and di2- where 
q 12 is defined by 

A(qi2, 612, 712) = A(ai2, 0 r 2 , 712) A t ( 0 3 , 63, 0 ), ( 37 ) 

or we could make the over-all Euler angles ay 2 , 0i2, and 712- It is equally easy to derive the 
matrix g^ for each case, but for illustration in this section, we consider the former case. 
This is also the situation considered in Refill, where it w r as called the Ijk coupling scheme. 
The later case wxmld be treated as in the next section. 

Thus w'e follow' the procedure of the previous section to obtain 
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for 0 = 1 and 2, and for 0 = 3, 
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where d[ 3 comprizes 7^, , and , wdiere BF is the body-fixed frame defined by 

the Euler angles 4>3, #3, and cv 12 . We then introduce the total angular momentum J BF to 


eliminate L BF : 


f BF _ _ fBF _ 7 B, 
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The functional form of jyf is given by rotating j , 2 by an- 
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to transform to differential expressions. Finally, we eliminate by using 
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In summary, we obtain 
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with g$ t given in Table III, and O t comprizing ^7, M 7, ^7, 37, ^7, 5777, J x > A ’ ana 

Jf F . It should be noted that when orthogonal coordinates are used, the rotation matrices 

A(0, A 2 , 712) will collapse to give the unit matrix when inserted into Eq. 8. Only for 

non-orthogonal coordinates do these give contributions to the functional form of the kinetic 

energy operator. 

Next let us re-address the issue of using 0 3, #3, and d, 2 rather than a, 2 , A25 an( l 712 
the overall Euler angles. Since one can show that (see Eq. 37 ) 


A(0, A 2 , 7 l 2 )A( 03 , # 3 ; “12) = A(Q,2, A21712), 


( 47 ) 
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if we introduce the rotation A(0, bj3 V2 , &712), then when 6 = 0, we recover the <p 3 , 0 3 , a l2 
case, while when 6 = 1, we recover the 0 * 12 , P 12 , 712 case. This is exactly analogous to the 
rotation A(O,a0 2 ,O) introduced in the previous section. Then we have 

A = £ 7 “) y A jk { 0, biin, l-in) £ sgo,, (48) 

"Aj J fc / 

where define the body-frame fXA As before, we have 


j* r = A(0 ,bp n ,bj n )J BF 


(49) 


but now 

— ^ > —ib sin bQri-J 1 ^ — ib cos b, 3 i 2 Jy J + - (^9) 

9/3i2 °P12 

and 

A -» -.wr + A. (51) 

c? 7 i 2 ‘ 0712 

We determined these last two relations by numerical experimentation. Thus we multiply 
g( 3 ) from the left by A(0,&& 2 , &7u), multiply the final three columns of the result from the 
right by A T (0, 6 ( 3 12 , 6712), then subtract ib sin 6712 times the fifth column from the seventh 
column, subtract ib cos 67 times the fifth column from the eighth column, and subtract ib 
times the sixth column from the ninth column to obtain the matrix g (3) in the expression 
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where O, comprizes m , m , efe' J * ' J v ’ 2 ‘ 

It should be noted that the rotation by 6/3 12 , 6y 12 has no effect on the kinetic energy for 


J = 0. 

Kinetic energy operators using the embedding defined by b have never been used before, 
and it would be very interesting to explore the benefits that arise from allowing 6 to be 
different from 0 or 1. One example that immediately comes to mind is this allows the 2 axis 
to be aligned with the symmetry axis for XH3 molecules at the equilibrium geometry. 

Although we do not do so here, it should be straightforward to intioduce Eckart rotations 
at this point. However, in contrast to the two vector case, it will not be possible to satisfy 
the Eckart relations by varying a single parameter. 
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VI. FOUR VECTORS 


When we have four vectors, one can use the results of a single vector and the results of 
three vectors, or one can use the results for two vectors twice. The former case is exactly 
analogous to the three vector case treated above. Following the pattern given before, one 
can easily derive the appropriate for either case. Since CH4 is a molecule of interest to 
us, we will only consider the case when we use two vectors tw ice. 

As in the three vector case, we need to change the twm vector notation. For the first two 
vector sub-part, w-e will call the angular momentum / 12 , the body-fixed frame defined for 
it bf 12, and the Euler angles ct vl , fin, and 7 12 . The quantities for the second two vector 
sub-part are defined analogously, except with the subscript 34. We will take ai 2 , /?i 2 , and 
7i 2 to define the body-fixed frame BF(= 6/12), so that we have for /5 = 1 and 2, 
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dx^ 


Y Aji{a 12) A2, 712) Y 9jJh 2^12’ 
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with () h , comprizing 5^-, 
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with () hA comprizing j b J 3 f , j%T , and jj§ 4 , with 

A ( A 34 , /? 3 4 , 734 ) = A(a 3 4,/?34 , 734)A T (Qfi 2 ,/?i 2 ,7l2)- 


(54) 


(55) 


We proceed as before by introducing the total angular momentum, which because jf/ and 
J BF are sign reversed w r hile is not, yeilds 


J,7 = ( 56 ) 

However j-ff is given by Eq. 42, so in the case of orthogonal coordinates wdiere the rotation 
A(ct34, ^34, 734) vanishes in the kinetic energy operator due to orthogonality, we obtain a 

r , , . r . d a d d d d aBF 

quite simple expression in terms 01 the 15 operators djl2 , dRi , 9Ra , dxi2 , dx3A < J x 34- 

jy34 ’ Jz3 4> , 3%\ fzM , 4? F , and J f F ' We Can US ° b ° th Eq ' 42 ^ Eq ' 43 t0 

convert j BF x etc to derivatives with respect to Q34, #34 , and 734, so there are only 12 unique 

operators. 

In Table IV we give the resulting matrix in symbolic form. There it is understood that 
the matrices g (21 ^ are the first three columns of g (2 \ g (2 ' 2) are the last three columns of g ! “\ 
and first three row’s are for j3 = 1, the next for /? = 2, etc. 
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However the kinetic energy operator derived in this way does not couple vibration and 
rotation in a symmetric manner. This is analogous to the triatomic case when one treats 
H 2 0 with a = 0. So we need to introduce the additional rotation A(ca 34 , c/3 34 , C734), which 
reverts to the above case when c = 0, makes the overall Euler angles a 34 , /? 34 , 7 34 when 
c = E and treats the two pairs of vectors more equivalently when c = 1/2. This rotation is 
included exactly as before, with the analog of relations 50 and 51 being 
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When one considers cases where c is nonzero, the resulting expressions can become quite 
complex. Since we are interested in including non-adiabatic corrections which in general 
requires numerical techniques, [17] we will take a numerical approach to evaluate matiix 
elements of the kinetic energy operator. In order to do this, one has to take care that the 
singularities are treated properly. We do this by writing 


g (4) = £ si 4 ’/., ( M ) 

S 

with f s equal to l/7?i, l/ 7 ? 2 , 1/7 1 3 , l/7? 4 , l/sinxi 2 i l/sinx 34 > an d 1/sin 8za- The g( ^ will 
always be finite. The we have 


G w = h 2 /2 £ f s fs' £ ( 61 ) 

55 ' jPS' 

and we use appropriate basis functions and quadratures which handle the singular factors. 


VII. DISCUSSION AND CONCLUSIONS 


We have discussed the procedure for producing kinetic energy operators for polyatomics 
using polyspherical coordinates and very general embeddings of body-fixed 2 axes. We do 
this by introducing the expression 

4- = £ 44 7 “) £ mO,- m 

OXi(3 j l 
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The matrices g can be very easily built up from the results from fewer numbeis of vectors, 
and the kinetic energy operator coefficients are simply given by dot products of columns of g. 
Noil-orthogonal coordinate systems pose no formal difficulty — however in practice, much 
more complicated expressions arise. We also show how to generalize the a parameterized 
axis embedding of Sutcliffe and Tennvson[7] to more than three atoms. 

One motivation for this work was the derivation of the expressions for the four vector 
case using non-orthogonal coordinates. One may ask why one would want to use non- 
orthogonal coordinates since the kinetic energy operator is so much simpler foi orthogonal 
coordinates? There are several reasons for this. One reason, put forward by Handy, [16] 
is that the potential matrix elements arc much simpler when one uses the physical bond- 
length-bond-angle coordinates. While less of an issue for triatomics, this becomes very 
important as the number of atoms increases, especially" if one wants to use a single potential 
energy surface for more than one isotopomer. The second reason is related to the treatment 
of Born-Oppenhcimer breakdown. W^e have showed the importance of including both first 
and second order corrections for H 2 O, [17] and since the second order correction gives rise 
to cross terms in the kinetic energy operator, even when orthogonal coordinates are used, 
the argument of simplicity no longer is valid. In addition, it is very useful for debugging 
programs to be able to compare the results obtained using different coordinates. Finally, it 
can be useful to turn to non-orthogonal coordinates to improve the convergence of the basis 
function expansion. [18] 

In addition to its utility for deriving expressions for the kinetic energy" operator, the ma- 
trices g are also required in generating non-adiabatic correction function originally 7 computed 
for cartesian nuclear position vectors. [19] 

In this work, we have not discussed what basis functions to use with these kinetic energy 
operators. Basically there are only two restrictions: first of all the functions must give finite 
matrix elements, and secondly they must adequately cover the Hilbert space for the problem 
of interest. For triatomics with a / 0, we used Jacobi polynomial basis functions to avoid 
the singularities in the matrix elements. [20] For the three and four vector case, it is not 
yet clear what changes one needs to make to ensure finite matrix elements when b or c are 
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non-zero. 


[1] M. Mladenovic, J. Chem. Phys. 112, 1070 (2000). 

[2] M. Mladenovic, J. Chem. Phys. 112, 1082 (2000). 

[3] F. Gatti, C. lung, M. Menou, Y. Justum, A. Nauts, and X. Chapuisat, J. Chem. Phys. 108, 
8804 (1998). 

[4] F. Gatti, C. lung, M. Menou, and X. Chapuisat, J. Chem. Phys. 108, 8821 (1998) 

[5] F. Gatti, C. Munoz, and C. lung, J. Chem. Phys. 114, 8275 2001). 

[6] M. Mladenovic, J. Chem. Phys. 113, 10524 (2000). 

[7] B. T. Sutcliffe and J. Tennyson, Int. J. Quant. Chem. 39, 183 (1991). 

[8] T. J. Lukka, J. Chem. Phys. 102, 3945 (1995). 

[9] V. Aquilanti and S* Cavalli, J. Chem. Phys. 85, 1355 (1986). 

[10] M. J. Bramley and T. Carrington, J. Chem. Phys. 99, 8519 (1993). 

[11] D. W. Schwenke, J. Phys. Chem. 100, 2867 (1996); 100, 18884 (1996). 

[12] A. R. Edmonds, “Angular Momentum in Quantum Mechanics”, third printing with correc- 
tions, Princeton University Press, Princeton, 1974. 

[13] R. T Pack, J. Chem. Phys. 60, 633 (1974). 

[14] O. Polyanski, poster at CCP6 workshop, Manchester, UK, Dec. 2001. 

[15] “Molecuar Vibrations”, E. B. Wilson, J. C. Decius, and P. C. Cross, Dover 1 ublications, Inc, 
New York, New York, 1955. 

[16] S. Carter and N. C. Handy, Mol. Phys. 47, 1445 (1982). 

[17] D. W. Schwenke, J. Phys. Chem. A, 105, 2352 (2001). 

[18] D. W. Schwenke, Chem. Phys. Lett. 189, 91 (1992). 

[19] D. W. Schwenke, to be published. 

[20] D. W. Schwenke, Comp. Phys. Commun. 70, 1 (1992). 


15 


TABLE I: Coordinate transformation matrix g (1) . 


A rbj j}f 

dR ^ x ^ y 

x 0 0 jji 

y o m 0 

z 1 0 0 
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TABLE II: Coordinate transformation matrix g' 2) . a = 1 - a. 


d 

ORi 

d 

dR‘2 

d_ 

Qx 

■h 

■h 

Jz 

— sin ax 

0 

— cos ax/ri 

0 

a cos ax/ iri 

0 

0 

0 

0 

- sinax/(*ri sinx) 

0 

— cosax/( ir i sinx) 

cos ax 

0 

— sin ax/i' l 

0 

asinax/'^i 

0 

0 

sin ax 

cosaxA’2 

0 

a cos ax/ ^2 

0 

0 

0 

0 

— sin ax/{ir 2 sinx) 

0 

cos ax/(ir 2 sin x) 

0 

cosax 

- sin ax/f '2 

0 

— a sin ax /i r 2 

0 
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TABLE III: Coordinate transformation matrix g^ 3 ^. a = 1 — a. 


d 

dRi 

d 

dRi 

d 

or 3 

d_ 

d 

dd 

dy 


Jy 

Jz 

— sin ax 

0 

0 

— cos ax/r i 

a cos ax 

—a cos ax sin 7 cot p/ 77 

0 

0 

a. cos ax sin 7 





x cos 7/77 




xcsc / 3 /zVi 

0 

0 

0 

0 

— sin ax sin 7 

-[sin ax cos 7 cot $ 

0 

0 

sinax cos 7 csc p 





xcsc xl r \ 

+ cos ax]csc x/ r t 



xcsc x/i r i 

cos ax 

0 

0 

— sinax/r’i 

a sin ax 

—a sin ax sin 7 cot p/r \ 

0 

0 

a sin ax sin 7 





X COS 7/77 




xcsc fi/ir 7 

0 

sin ax 

0 

cos axh'2 

a cos ax 

—a cos ax sin 7 cot p/r 2 

0 

0 

a cos ax sin 7 





x cos 7/7*2 




xcsc P / il '2 

0 

0 

0 

0 

— sin ax sin 7 

— [sin ax cos 7 cot 8 

0 

0 

sinaxcos^csc P 





xcsc x/ r 2 

— cos ax]csc xh 2 



xcsc x/^ r 2 

0 

cosax 

0 

— sinax/ r 2 

-a sin ax 

a sin ax sin 7 cot Pjr 2 

0 

0 

—a sin ax sin 7 





x cos 7/7*2 




xcsc P / 77*2 

0 

0 

0 

0 

-1/^3 

0 

0 

l/*>3 

0 

0 

0 

0 

0 

0 

CSC /?/t*3 

— l/l >'3 

0 

— cot /3/ir3 

0 

0 

1 

0 

0 

0 

0 

0 

0 
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TABLE IV: Coordinate transformation ma :rix g (4) ■ 

a a a a a a a a a j j j 

dR\ dR'2 dxi2 dR$ d#4 #X34 da 3 4 d0z4 ^73 4 x y 

^TT 0 

0 A r g( 21 ) A T g( 2 - 2 )R 2 0 

cos 0:34 cot /?34 sin 0:34 — cos a^csc$-n 

R! = sin 0:34 COt ^34 —COS 0:34 — sin d34CSC/?34 

-10 0 

— COS 734 CSC/334 siri 734 COS 734 cot /?34 

R '2 = sin 734 csc/ 934 cos 734 — sin 734 cot /I34 

0 0 1 

A 7 ’ = A T (d 34 ,^34,734) 0 

0 A 7 (0:34, ^34, 734) 
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